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ABSTRACT 

When applied to deep photometric catalogs, the two-point angular correlation function, uj{9), is a 
sensitive probe of the evolution of galaxy clustering properties. Here we present measurements of 
uj{0) as a function of Iab magnitude and {R — I) color to a depth of Iab = 24 on scales 6 = 7" — 3' , 
using a sample of ^ 350, 000 galaxies covering 5 degrees^ in total over 5 separate fields. Using redshifts 
of 2954 galaxies in early DEEP2 Galaxy Redshift Survey data, we construct robust galaxy redshift 
distributions as a function of /ab and Rab magnitude and {R—I) color for galaxies between < z < 2. 
We constrain models of the redshift evolution of galaxy clustering and find that significant growth 
of clustering has occurred from z > 1 to z = 0. A model in which the comoving scale- length, xq, 
ev olves line arly with redshift, xo{z) = 2:0(0) (1 — Bz), fits the data better than the e model proposed 
bv lGroth~fc Peebles (JL977). The clustering properties depend strongly on observed {R — I) color, 
with both the reddest and bluest galaxies exhibiting large clustering amplitudes and steeper slopes. 
Different observed (i? — /) color ranges are sensitive to very disparate redshift regimes. Red galaxies 
with (R — I) ~ 1.5 lie in a narrow redshift range centered at z ~ 0.85 and have a comoving scale 
length of clustering at z = 0.85 of xq = 5.0±0.3 Mpc. These galaxies have early-type spectra and 
are likely progenitors of massive local ellipticals. The bluest galaxies with (R — I) ^ appear to be a 
mix of star-forming galaxies, both relatively local (z ^ 0.3 — 0.6) dwarfs and bright z > 1.4 galaxies, 
and broad-line AGN. We find that local blue dwarfs are relatively unclustered, with xq = 1.6 ± 0.2 
Mpc. The z > 1.4 blue galaxies have a larger clustering scale-length, 0:0 >5 Mpc. 

Subject headings: galaxies: statistics — galaxies: evolution — cosmology: large-scale structure of 
universe 
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1. INTRODUCTION 

Studies of galaxy clustering can constrain cosmological 
parameters as well as galaxy formation models, since the 
observed galaxy clustering strength depends on both the 
underlying dark matter distribution and the efficiency of 
galaxy formation inside dark matter halos. While the- 
ory and simulations both provide detailed predictions 
for the dark matter distribution in various cosmological 
models, the processes of galaxy formation and evolution 
are less well-understood. The redshift dependence of the 
galaxy clustering strength, in particular, can be used to 
constrain galaxy evolution theories if the cosmology is 
known. 

Galaxy clustering is ideally studied in three dimen- 
sions, where the redshift of each galaxy is known. How- 
ever, cosmic variance dictates the need for very large 
data samples, probing a large volume of space and sam- 
pling multiple independent fields, in order to obtain sta- 
tistically robust results. As photometric data exist over 
larger areas and to fainter magnitude limits than are cur- 
rently available for spectroscopic samples, galaxy cluster- 
ing has long been measured from two-dimensional pho- 
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tometry alone. 

The angular two-point correlation function, uj(9), mea- 
sures the excess probability above random of finding a 
galaxy at a specific angle 9 from another galaxy. This 
quantity is derived from counts of galaxy pairs as a func- 
tion of their separation, and while relatively straightfor- 
ward to measure, its interpretation is not trivial. From 
the projected angular two-point correlation function, 
uj{0), one can infer the three-dimensional spatial two- 
point correlation function, ^(r), if the redshift distribu- 
tion of the sources in each sample is well-known and the 
redshift dependence of the clustering strength is known 
or assumed. The two-point correlation function, ^(r), is 
usually fit as a power law, ^(r) = if/'''o)~'^, where ro is 
the scale-length of the galaxy clustering. However, the 
assumed galaxy redshift distribution {dN/dz) has varied 
widely in different studies, such that similar observed an- 
gular clustering results have led to widely different con- 
clusions. A further complication is that each sample usu- 
ally spans a large range of redshifts and is magnitude- 
limited, such that the mean intrinsic luminosity of the 
galaxies is changing with redshift within a sample. As 
the clustering properties of galaxies are expected to de- 
pend on luminosity and have been shown to do so locally 
|(Nprbcrg et al. 2001), this can hinder interpretation of 
the evolution of clustering measured in iv{9) studies. 

Many of the first useful measurements of large-scale 
structure were studies of angular clustering. One of 
the earliest measurements of the uj{9) statistic was the 
pioneering work of Peebles (,1975) using photographic 
plates from the Lick survey. They found uj{9) to be 
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well fit by a power law, uj{9)= A;^9^ , with a slope of 
S = —0.77 ± 0.06. Later studies using CCDs were able 
to reach deeper magnitude limits and found that fainter 
galaxies had a lower clustering amplitude. However, 
most of these surveys covered only one or two fields, 
typically with a field of view of ~ 50 arcmin^ due to 
the small sizes of C CDs, and were strongly affected by 
cosmic varian ce (e.g., lEfstathiou et alJl99 lURoch e et all 
El93; Braine rd. Small. & Moul dl 119951: iH udon fc Lillvl 
IT996t IWoods fc Fahlman 19971E iMcCracken et al. 2000). 

More recently, wide-field cameras on intermediate- 
class telescopes have led to deep wide-field photo- 
metric surveys (e.g., LN euschaeferJ^_ffiindhQrst! '1995!; 
[ Roche fc Ealesll999blCab anac. d e Lapparent. fc Hickson 
boOOt IMcCracken et al.ll2001t (Wilson ^OOS") . A ground- 
breaking study was conducted by Postman et al. (1998), 
who used a 4-m telescope with a 16 arcmin field of view to 
cover a much wider area by tiling together several point- 
ings, surveying a contiguous 4 x 4 degree^ field to a depth 
of Jab = 24, reaching to z ~ 1. IMcCracken et al.l lj200lD 
present a survey of 1 degree^ to a depth of Jab = 25.5 
over four in dependent fiel ds using the UH8K camera on 
CFHT, and lWilsoul l|200l use the same camera to image 
1.5 degrees^ over 3 separate fields to a depth of / = 24.0. 

Other recent studies employ accurate photometric red- 
shifts to measure w{9) for galaxy sampl es of a spe- 
cific ty pe over well-defined redshift ranges. iBrown et al.l 
([2003(1 investigate the clustering of red galaxies in 
three redshift bins between 0.3 < z < 0.9 in a 1.2 
degree^ field from the NOAO Deep Wide-Field Sur- 
vey pannuzi fc Dev 1999), while Budavari ct al. (2003) 
use data from t he Sloan Digital Sky Survey (SDSS, 
lYork et all l[2000D ) to construct a volume- limited sample 
with more than 2 million galaxies between 0.1 < z < 0.3 
with which to study galaxy clustering as a function of 
color and spectral type. 

In this paper we present a survey undertaken with the 
CFH12K camera covering a total of 5 degrees^ over 5 
separate fields to a depth of Iab — 24.0. There is a wide 
diversity of adopted redshift distributions for magnitude- 
limited samples in the literature, which can lead to dif- 
ferent interpretations of similar results. Here we use 
much improved dN/dz distributions, relying on recent 
data from the DEEP2 Galaxy Redshift Survey. Unlike 
previous studies, we require no spectral modelling in fit- 
ting for these dN/dz distributions and no assumptions 
are made about the evolution of the galaxy population. 

An outline of the paper is as follows. In Sections 2 
and 3 we describe the observations and data reduction 
process. As a large subset of this data was used for the 
DEEP2 Galaxy Redshift Survey photometry, we provide 
a detailed account of the data reduction pipeline. In 
section 4 we present galaxy counts in three photometric 
bands, while in Section 5 we measure uj{9) as a function 
of Iab magnitude, describe the dN/dz distributions we 
use and discuss constraints on the evolution of galaxy 
clustering. Section 6 parallels Section 5 but focuses on 
clustering as a function of (i? — /) color. In this section 
we discuss in particular the strong clustering seen for the 
reddest and bluest galaxies and discuss the most likely 
interpretations of this effect. We conclude in Section 7. 



Our data were taken using the CFH12K camera on the 
3.6-m Canada-France-Hawaii telescope. The CFH12K 
camera has a 12,288 x 8,192 pixel CCD mosaic array and 
a plate-scale of 0.207" per pixel, providing a field of view 
of 0.70 X 0.47 degrees. We observed five separate fields 
on the sky, with one to five distinct CFHT12k pointings 
per field, where the area in one pointing equals the field of 
view of the camera. Our data were taken during several 
observing runs from September 1999 to October 2000 
and with queue-scheduled observing from April to July 
2001; observing details are given in Tabled 

We integrated on each pointing for ~1 hour in B and R 
and ~ 2 hours in /, broken down to individual exposures 
of 600 seconds. We dithered the telescope slightly be- 
tween exposures to fill in gaps between the CCDs and 
to cover areas with low charge transfer to yield well- 
sampled combined images. We also took four additional 
120-second /-band astrometry images in each pointing, 
with offsets of 3.5' in RA and 7' in Dec. between images. 
The seeing varied from 0.69" to 0.95", with final seeing 
estimates for each pointing listed in Tabled 

3. DATA REDUCTION 

3.1. Initial Processing 

The data were reduced using the C-based IMCAT soft- 
ware package developed by Nick Kaiser^. After bias- 
subtracting each exposure, we constructed a super-flat 
for each passband from the median of the science im- 
ages, which were individually normalized by the mode. 
The super-flats were adjusted by applying a multiplica- 
tive correction to each chip (aside from one chip arbitrar- 
ily chosen as a reference) such that when divided into the 
science data frames they resulted in a near constant sky 
level across the entire focal plane. We then created mask 
images which masked out saturated stars, regions with 
bad charge transfer and other CCD defects. After di- 
viding each exposure by the appropriate super-flat we 
applied the mask images. A model for the sky level in 
each exposure was then obtained by estimating the local 
mode using a set of 8x16 sub-images of size 256x256 
pixels and smoothly interpolating between these points 
using triangular tessellation (very discrepant values, e.g. 
around bright stars, were rejected). We then created ini- 
tial catalogs of objects from each of the sky-subtracted 
images, computing magnitudes using a 20-pixel radius 
aperture, which is large enough to be unaffected by vari- 
ations in the seeing, but not so large that systematic 
errors in flux from uncertainty in the background sky 
level dominate. From these catalogs we selected stars in 
each exposure, for both the astrometric and science im- 
ages, selecting compact, bright objects. The star catalogs 
from the astrometry images were then used to determine 
astrometric solutions for our dataset. 

3.2. Astrometry 

We used the US Naval Observatory catalog (USNOA 
2.0) as an absolute astrometric reference system. In each 
pointing, we used stars which were detected in both the 
USNOA catalog and in at least one of our four astrome- 
try images to determine a set of true star positions from 
which we fit (rejecting outliers) cubic polynomial distor- 
tion parameters for each of the 4 exposures x 12 chips. 



2. OBSERVATIONS 
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The resulting rms errors are approximately 0".006, as 
shown in Figure ^ where we plot the residuals in pixel 
coordinates of the stellar positions from the best-fit so- 
lutions. These are our relative astrometric errors; our 
absolute errors are defined by the accuracy of the US- 
NO A catalog, which is roughly 0.5". 

We then used these field distortion models to compute 
the amplification associated with the distortion as a func- 
tion of position on the focal plane. This can in turn be 
used as a check of our astrometry solutions; we expect 
the amplification to be a simple quadratic function of 
radius, such that gross deviations from this form would 
most likely indicate systematic errors (feeding in from 
the large-scale systematic errors in the USNOA catalog) . 
The field distortion amplifications were also needed to 
correct the stellar magnitudes. We found that the field 
distortion conformed well to the expected circularly sym- 
metric quadratic behavior, and we conclude that system- 
atic errors in our astrometry are very small. After solv- 
ing for the astrometric solution between the astrometry 
images and the USNOA catalog, we fit for astrometric 
solutions between the science exposures for each point- 
ing and the astrometric exposures, again using catalogs 
of stars selected by size and magnitude. 

3.3. Initial Photometry 

Each science exposure yielded a catalog of ^ 800 — 1000 
stars. The stellar magnitudes were corrected for field- 
distortion amplification, and the catalogs of each pair of 
exposures were merged using the sky coordinates pro- 
vided by the astrometric solution. This provided a very 
large set of pairs of measured stellar magnitudes, along 
with the chip numbers, detector coordinates and expo- 
sure numbers. We then fit, by least squares minimiza- 
tion, a model in which the magnitude of an object mea- 
sured on chip c in exposure e is mce = m + drric -\- drue- 
This method simultaneously measures corrections for ex- 
tinction terms in individual exposures (dnie) as well as 
chip-to-chip zeropoint offsets {drric). The extinction val- 
ues fit were generally very small, typically 0.005 mag- 
nitudes, indicating that conditions were accurately pho- 
tometric, but much larger values were found for the chip 
correction factors dmc, ranging up to ^ 0.05mag. This 
is not surprising, as the CCDs in the CFHT12K are an 
inhomogeneous collection, with substantial variations in 
the quantum efficiency. The resulting error at bright 
magnitudes (16.5 < m < 18.0) is ((mi - ma)^)^/^ « 0.02 
mag, as shown in Figure |21 At fainter magnitudes, the 
scatter is larger due to measurement error. These results 
show that with this simple correction the magnitudes for 
stars are highly reproducible cxposure-to-exposure with 
an error - 0.02/\/2. 

Finally, we combined the individual science exposures 
for a single pointing, deleting cosmic rays, correcting for 
the field distortion, and creating a set of 9 x 13 sub- images 
of size 1024 x 1024 pixels ('quilts'). For each quilt, we 
first determined which input images contributed; then 
the appropriate sub-images were generated and combined 
into a set of image planes. These were median-combined, 
with rejection of discrepant pixels. 

3.4. Catalog Construction 

Object catalogs were constructed for each quilt image 
using the faint object detection algorithm described in 



Kaiser, Squires, & Broadhurst (1995). We smoothed 
the i?-band images (the deepest obtained) with a se- 
quence of 'Mexican hat' filters of progressively larger 
size, tracked the peak trajectories through the three- 
dimensional space of X, y, and Gaussian radius, and then 
defined an object to be the peak of the significance (the 
height of the peak divided by the rms fluctuation for that 
smoothing scale) along the trajectory. We used a filter 
which is a normalized Gaussian of scale rg minus a nor- 
malized concentric Gaussian of scale 2rg. The filter scale 
was varied from 0.5 to 20.0 pixels with equally spaced log- 
arithmic intervals Ar^ — 0.2rg; the best-fit value of rg 
provides an excellent estimate of the Gaussian radius of 
an object on the sky, with smaller fractional scatter (for 
stars) than either the half-light or Petrosian radii mea- 
sured by IMC AT. We therefore use Vg as our measure 
of object size throughout the remainder of this paper. 
Note that because of the hmits of the grid, objects with 
rg > 20 pixels will be poorly measured; however, this is 
of little consequence in the magnitude ranges considered 
here. All objects with significance ly > A were output to 
an initial object catalog; we then perform photometry for 
each object in four different circular apertures (of radius 
3rg, 6rg, 1", and 2") in each band. 

Unfortunately, the algorithm used to find objects can 
produce false detections due to instrumental features 
such as diffraction spikes, bleed trails, regions of scat- 
tered light, etc. To remedy this, we developed separate 
IDL software to identify the problematic pixels on the 
projected, combined CFHT12k images; this is done sep- 
arately in B, R, and /. Our algorithm begins by flagging 
all saturated pixels, along with the central pixels of all 
partially-saturated sources identifiable in the initial ob- 
ject catalog; this limit is found empirically in each point- 
ing as the magnitude where the measured sizes of stars 
are no longer independent of their flux. Then, all pixels 
less than 3 pixels from a flagged point in both X and Y for 
which the flux (after smoothing with a 3-pixel boxcar) is 
more than 3.25 times higher than the sky noise (also after 
smoothing, from the IMCAT sky inverse variance map) 
are themselves flagged as bad. This procedure is iter- 
ated until there are no pixels with significant flux within 
3 pixels of a flagged pixel. In what follows, all objects 
whose centers are within 10 pixels of any flagged pixel (or 
pixel with no data, e.g. due to edges or bad pixels in the 
detector) in any band are rejected from analysis; this ra- 
dius corresponds to ^ 4rg for a typical object, and thus 
is larger than the flnal apertures used. The bad pixel 
masks produced by this procedure are also employed in 
masking the random catalogs used for correlation analy- 
ses (see Section 5.1). 

We found by visual inspection of the resulting masked 
catalog that occasionally there were spurious or incor- 
rect entries, either due to overlapping objects, where 
substructure within a large, extended object is detected 
as a small separate object within its bounds, or "over- 
merged" pairs or triplets of galaxies which were identifled 
incorrectly as one extended object. To rectify these sit- 
uations, we applied two algorithms, one to delete over- 
lapping objects and one to split merged objects. For 
the overlapping objects, we first define a 'small' radius, 
rg = 0.85rseeing, and then flag as spurious any objects 
with rg < Ts and which are inside 3rg of a bright, ex- 
tended object or that have rg < 1.5 * and are within 
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3rg of a bright, extended source and have peak counts 
more than 2 sky sigma below the peak of the nearby 
background (i.e. he on a very strong gradient back- 
ground). For all tests of object overlap, we use ellipsoids 
with rg/ivg as the semimajor axis, which matched vi- 
sual morphologies well. As a part of this routine, we also 
eliminate from the catalog identical overlapping objects 
identified on multiple quilts. 

For the over-merged objects, we first fiag all low 
surface-brightness (LSB), elongated objects, with b/a < 
0.5, as many of these are multiples which have been 
incorrectly designated as a single object. Around each 
flagged object we smooth the image quilts with a Gaus- 
sian kernel of cr = 0.67aseeing and find all peaks within 
a 2-pixel radius circle which are greater than five times 
the sigma of the smoothed sky (i.e., the highest surface 
brightness regions present). We then centroid on these 
and measure their FWHM on the image quilt. We evalu- 
ate the significance of the original object by determining 
whether, after masking the candidate sub-objects, a fit 
to an elliptical-Gaussian model has an amplitude greater 
than 5ct above (where here a is the uncertainty in this 
Gaussian- weighted amplitude). If not, we give the origi- 
nal object a quality flag such that it will be rejected from 
further use. If the signiflcance of detection for the origi- 
nal object is > 30, then we consider all the high surface 
brightness (HSB) regions to be merely bright sub-regions 
of an original object, and reject them. Otherwise, for all 
the HSB regions, we evaluate: 1) whether they are resolv- 
able from other HSB regions identified in this process - in 
case of such conflicts, we keep only the brightest peak; 2) 
if the HSB object is indistinguishable from a known ob- 
ject in the catalog - i.e. we require its distance from each 
known catalog object be: a) > 1.25aseeing from the cat- 
alog object, such that separation of the peak is resolved, 
or b) 0.67 times the catalog object profile a from central 
object, so that we are not likely to be seeing substruc- 
ture in the peak; and 3) whether the identified object 
in fact appears diffuse - i.e., has FWHM> 2.5 times the 
seeing FWHM, and not identified with a catalog object. 
This throws out false diffuse detections due to low-level 
flat-field or sky-subtraction issues that are not the over- 
merged objects that we are trying to eliminate. All of 
the criteria and limits used were tuned empirically based 
on the actual over-merging of objects seen. We found 
that ~ 3% of the objects were "over-merged" and were 
thus split into separate objects, on which we performed 
photometry in all three bands. As colors are computed 
in 1" apertures, biases in the colors due to blending of 
sources are minimized. 

3.5. Star- Galaxy Separation 

The magnitude range of our data spans a regime where 
stars are a non-negligible contaminant for a galaxy cat- 
alog. Unfortunately, size and magnitude alone are not 
always sufficient to distinguish stars from galaxies, es- 
pecially at fainter magnitudes where small intrinsic sizes 
and uncertainties in size measurements cause some galax- 
ies to have measured sizes comparable to those of stars. 
We have therefore developed a set of algorithms which 
distinguish stars from galaxies in a probabilistic fashion, 
using their sizes, magnitudes, and {B — R), {R — I) col- 
ors. We exploit the fact that stars and compact galaxies 
have very different magnitude distributions, with stars 



dominating at bright magnitudes and galaxies at faint 
magnitudes, and that stars fall within a tight locus in 
color-color space. Objects which are clearly extended, 
with rg > 0.05 above the upper limit of the locus of stars 
in the size-magnitude plane, are deemed galaxies, while 
compact objects are given a probability, in a Bayesian 
sense, of being a galaxy given their magnitude m and 
color vector c: 

p i{m c) = P{c\gal)p{gal\m) 

^^"^ ' p{c\gal)p{gal\m) + p{c\star)p{star\m) 

The magnitude distribution of compact objects is es- 
sentially linear at the bright end (i? < 21) and rises ex- 
ponentially at fainter magnitudes. We assume, quite rea- 
sonably, that the linear component is due exclusively to 
stars, while galaxies cause the exponential increase at 
the faint end. In each pointing, we perform a non-linear 
least-squares fit to the magnitude distribution function 
n{m) ~ a^e"'^"^ + a2 + a3m, where m = mii — 22. We then 
define the probability of a compact object being a galaxy, 
based on its magnitude, as p{gal\m) = (aoe"i™)/n(m), 
with p{star\m) = 1 — p{gal\m). In practice we limit the 
range for these priors, confining them to the range 5%- 
95%, in order to restrain the influence of any errors in 
our assumptions. It then remains to determine the color 
distributions of stars and galaxies. 

It is first necessary to bring all pointings to a common 
color system. We do an initial, rough calibration using 
observed number counts; in the absence of cosmic vari- 
ance, the number of galaxies within a given magnitude 
range should be the same in each pointing and bandpass. 
We can refine the zeropoints of our color system by forc- 
ing the locus of stars, chosen for their compact size and 
bright magnitudes, for the first iteration of this proce- 
dure, in color-color space to match amongst all point- 
ings. This locus is well-approximated in our bands by 
two intersecting line segments of fixed slope; therefore, 
by fitting the color distribution of stars in each pointing 
to that functional form, via least-squares minimization, 
we may derive the offset in both {B — R) and {R — I) col- 
ors required to make the stellar locus match our reference 
pointing. We chose pointing 2 in field 2330 (see Table ^ 
to be our reference pointing, for reasons described be- 
low. This procedure indicates that the color zeropoints 
derived from number counts are consistent from point- 
ing to pointing to within 0.03 mag rms in (i? — /) and 
0.05 mag rms in {B — R). By using the stellar locus, 
although we do not necessarily know the correct abso- 
lute zeropoint in any pointing in any band, we can bring 
all pointings to a common color system in a robust way. 
We iteratively repeat this stellar locus calibration as the 
color distributions of stars and galaxies are refined, using 
the probabilistic star-galaxy separation defined above af- 
ter the first color distributions are constructed; excellent 
convergence occurs within two iterations. In addition, we 
use the variations in the stellar locus within subregions 
of each individual pointing to derive a smooth correction 
for residual scattered light and/or flat field errors, ap- 
plied using a minimum-curvature surface interpolation. 

After bringing all pointings to a common color sys- 
tem, we select objects within a pointing which arc: 1) 
almost certainly stars (i.e., that have rg within the stel- 
lar range, and either 19 < m/j < 21.25, or mji < 24, 
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{R — I) > 1.45, and color within ±0.2 mag of the ex- 
pected stellar locus - this latter requirement is neces- 
sary to include M stars), 2) bright compact galaxies 
(20.5 < mR < 23.5, rg larger than the stellar range 

but less than ^{l.bpix)'^ + r'^eeing)^ '^^ fa-i^t compact 

galaxies (m_R > 23.5 and Vg as above). Plots of the po- 
sitions for compact galaxies and stars in (B — R) and 
(i? — /) space for data in one pointing are shown in Fig- 
ure |3| 

We then create a two-dimensional color-color his- 
togram for each of these samples, which we smooth by 
a two-dimensional Gaussian kernel with a cr = 0.05 
mag to take into account the effects of magnitude er- 
rors in applying these distributions to catalog objects. 
We delete outliers, redo the smoothing, and normalize 
the distributions to have a maximum probability of one. 
This results in a probability distribution for the color 
of stars directly, though for galaxies we have to take 
into account the magnitude-dependence of their distri- 
bution in color-color space, so as not to be dominated by 
the more numerous, faint objects. We therefore define 
p{c\gal) = w(m}i)p{c\hright) -I- (1 — w{mfj))p{c\faint), 
where bright and faint denote the two galaxy samples 
we defined above and w — 1 for mji < 23, w = uir — 23 
for 23 < mu < 24 and w = for uir > 24. We then use 
equation2]to determine the probability for each compact 
object that it is a galaxy. 

We emphasize that we have computed color-color dis- 
tributions for the bright and faint compact objects sep- 
arately; we are not assuming that galaxy colors and 
magnitudes are completely independent. It is primar- 
ily important in our Bayesian method that the color 
distribution of galaxies - our prior - fills a different 
space than the stellar color distribution. And as galax- 
ies of a given redshift and restframe color have substan- 
tial breadth to their luminosity function, the covariance 
between color and magnitude is dilute. A full treat- 
ment of the three-dimensional distribution of magnitude 
and color would greatly complicate implementation and 
would yield noisy results at bright magnitudes where we 
have fewer galaxies to build a distribution from, but also 
where star-galaxy separation is most important. The re- 
sulting strong separation between stars and galaxies that 
we find using our probabilistic method shows that we are 
not sensitive to such details. 

We find that 68% of the objects are extended, while 
2% do not have sufficient information to determine Pgah 
either because they lack data in one or more passbands 
or are brighter than the saturation limit. Figure ^ shows 
the Pgai distribution for the compact objects, which is 
extremely bi-modal, with 25% of the objects having 0% 
probability of being a galaxy and 38% having 100% prob- 
ability of being a galaxy. The few objects with intermedi- 
ate values of pgai tend to lie near these extremes, and so 
in the analyses of angular correlations that follow we in- 
clude only those objects with a > 20% probability of be- 
ing a galaxy, though it makes very little difference what 
exact cutoff is used. Effects on the angular clustering 
measurements due to stellar contamination are discussed 
in section 5.2. 

3.6. Final Photometry 



A large portion of the photometry prese nted here is 
used b y the DEEP2 Galaxy Redshift S^irvev lPavis et alJ 
()2002j) , and as a courtesy we define the photometric cal- 
ibrations such that they are optimal for that survey. We 
therefore define a set of magnitudes in the AB system, 
with a total, i?AB-band magnitude and {B — R), {R — I) 
colors for each object. We adopt as our measure of to- 
tal magnitude the R magnitude measured in a circular 
aperture of radius 3rg, unless 3rg is less than 1" (due to 
noise in its measurement), in which case we use the mag- 
nitude within a 1" radius. We found that the magnitude 
within 6rg as measured by IMCAT is often contaminated 
by nearby objects included within the aperture or prob- 
lematic as the aperture goes beyond the buffer region 
of each quilt for objects near their edges, while the 3rg 
magnitude proved more robust. 

Conventionally, we would then use the B and / magni- 
tudes measured within this same aperture to define our 
system; however, due to differences in seeing/aperture 
corrections between bands and amongst the CFHT12k 
pointings, it would then be impossible to make the col- 
ors of extended objects (i.e. galaxies) agree pointing-to- 
pointing even if the colors of the stellar locus match; this 
effect can reach ^ 0.05 mag in each band. This effect 
is much smaller for colors derived from the IMCAT 1" 
aperture magnitudes, < 0.02 — 0.03 mag in both {B — R) 
and {R — I), based on cross-correlating the color locus of 
extended objects from pointing to pointing. Therefore, 
we adopt as the basis of our magnitude system the quan- 
tities we can measure and calibrate most robustly: the 
"total" Rab magnitude, measured in the 3rg aperture, 
and the two colors {B — R) and {R — I), each measured 
within 1" apertures. In the remainder of this paper, B 
will be used to represent Rtotai + {B — R)^ii, and / will 
represent Rtotai — [R— I)^"- Note that these should be 
distinct from conventional magnitudes only in the case 
of resolved color gradients within the observed galaxies. 

We calibrate our i?AB-band magnitude and {B — R), 
{R — I) colors through two separate processes. First, 
we calibrate B, _R, and / to an AB magnitude system 
in the native CFIIT12k passbands using SDSS photome- 
try for stars in our reference pointing (pointing 2 in field 
2330 [see Table ^). This pointing was chosen to define 
our color system as it yielded the cleanest comparison to 
SDSS, due to the quality both of the CFHT12k observ- 
ing conditions and the SDSS photometry (it hes on the 
celestial equator and has been observed by SDSS repeat- 
edly). The accuracy of this calibration is ~ 0.02 — 0.03 
mag (worst in B), dominated by errors in the SDSS cal- 
ibration. Through iteratively fitting the stellar locus, 
bringing all pointings to a common color system, deter- 
mining the overall color distributions of stars and galax- 
ies, and reapplying our star-galaxy separation algorithm, 
we not only produce a stable star-galaxy separation but 
also bring all pointings to the same {B — R)is, {R — I)is 
color system to within ~ 0.01 mag in each color; the ze- 
ropoints of that color system are then set by the SDSS 
calibration of our reference pointing. After our correction 
for spatial variations in the color system using variations 
in the stellar locus, residual variations in the {B — R) and 
(i? — /) zeropoints within each pointing are estimated 
to be 0.04 and 0.03 mag, respectively; these variations 
should primarily be relevant on scales of ~ 4' — 8'. 
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Second, after this calibration of the color system has 
been performed, we calibrate the i?-band zeropoint in 
each pointing. In the four pointings for which no infor- 
mation from SDSS is available (and which are not used 
in the DEEP2 Galaxy Redshift Survey), we determine 
zeropoints by forcing the number counts to match those 
from our other fields; based on pointings with SDSS cal- 
ibration, we find that this yields a zeropoint with errors 
of 0.10 mag, due to cosmic variance (errors from this 
technique in {B — R) or {R — I) are smaller as the cos- 
mic variance is strongly covariant between bands). We 
have tested several of the below analyses with and with- 
out the inclusion of these less-calibrated pointings, and 
all differences were well within the quoted errors. In 
the fields with SDSS calibration, however, by comparing 
B, R, and / to the SDSS magnitudes for stars, we get 
three independent estimates of the R zeropoint (as we 
have tied the three magnitudes together in our photom- 
etry by making the stellar loci match). In general, the 
zeropoints estimated from R and / agree to within 0.02 
mag rms, while these agreed with the zeropoint from B 
to ^ 0.04; this matches well the expected variation in 
SDSS photometric zeropoints. We use the average of the 
zeropoint offsets from R and / to correct our photometry 
wherever SDSS photometry is available. Based on com- 
parisons to SDSS, residual variations of the R zeropoint 
within each pointing, which we cannot correct for with 
either the stellar locus or number counts, appear to be 
< 0.02 mag rms. The final uncertainty in the absolute 
calibration includes both a random error due to the finite 
number of stars with well-measured magnitudes in both 
datasets and a systematic error due to uncertainties in 
the AB magnitude calibration of the SDSS photometry. 

4. GALAXY COUNTS 

We first present galaxy number counts for our data in 
order to test the uniformity and depth of the sample. 
Figure El shows histograms of Bab, ^ab and Iab galaxy 
counts for all 15 pointings, in units of log number of 
counts per degree^ per 0.5 magnitude. At the bright 
end we see variations in the counts due to cosmic 
variance, while at the faint end the varying depth of the 
data in each pointing leads to differences in the counts 
past our completeness limit. Our data is complete to 
~ 25.25 in Bab, ~ 24.75 in i?AB, and - 24.25 in Iab 
in all pointings except field 1416, pointing 2, which 
is our shallowest pointing and the low data point at 
each of the completeness limits in Figure El We fit 
for the slope in the galaxy counts, 6 = d logN/dm, in 
each band over the magnitude range in which the data 
is consistent from pointing to pointing. The B-band 
galaxy counts have a slope oi 5 = 0.49 ± 0.01 in the 
magnitude range Bab = 21.25 — 24.75 mag. This fit is 
shown in the bott om of Figure O as the dotted l i ne and 
agrees well with IBrunner. Cqn noUv. fc Szalav* flDO!/) 
{S = 0.51, Bab ~ 20 - 24), [McCracken ct al. (2001) 
(6 = .47 ± 0.02, Bab = 20 - 24) and iMcCracken et all 
112001 {5 = 0.45 ± 0.01, Bab = 20 - 24). Our i?AB 
galaxy counts have a slope oi S — 0.37 ± 0.01 in the 
magn itude range i? a b = 2 0.25 — 24.25 mag, in agreement 
with Ho gg et aP lll997l ). IBru nner. C onnolly. & Szalav 
119991. iMetcalfe et all l|200H) and IMcCracken et all 
(200l who find slopes of (5 = 0.33,0.34,0.37, and 
0.37 ± 0.03, respectively. Our Iab counts have a 



slope of ^ = 0.33 ± 0.01 in the magnitude range 
^ab = 19.75 — 23.75 mag, which also agrees well 
with other published results in similar magnitude 
ranges, all of whic h report a slo pe be twe en 0.31 — 0.35 
iBnmner. ConnoU v. fc Szalavl II1999D. | Arnouts et al 
Il999), Metca lfe et alJ ll200lD. IMcCracken et al.. (,2001l 
IWllsonl (120031) .' IMcCracken et all i200.f T 
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GALAXY CLUSTERING AS A FUNCTION OF 
MAGNITUDE 



In this section we present measurements of lu{9) as a 
function of Iab magnitude and use redshift distributions 
of galaxies from the early DEEP2 Galaxy Redshift Sur- 
vey data to constrain models of clustering evolution. Sec- 
tion 6 presents a parallel study of galaxy clustering as a 
function of {R — I) color. 

5.1. The Angular Two-point Correlation Function 

The projected angular two-point correlation function 
uj{9) is a measure of the probability above Poisson of 
finding two galaxies with a separation 9, defined as 



dP = N[l + uj{9)]dn 



(2) 



where TV is the mean number of galaxies per steradian 
and dfi is the solid angle of a second galaxy at a sep- 
aration 9 from a randomly chosen galaxy. To measure 
to (9) one must first construct a catalog with a random 
spatial distribution and uniform density of points with 
the same selection criteria as the data, to serve as an un- 
clustered distribution with which to compare the data. 
For each of our 15 pointings we created a random cata- 
log with 200,000 points (17 times larger than the largest 
data sample we use here) with the same sky coverage as 
our data, for which we then applied the same masking 
that was used to mask bad regions of the data due to 

saturated stars and CCD def ects. 

We measure u!{9) using the lLandv fc Szalavl ljl993|) es- 
timator, 



u;(9) = — 
^ ' RR 



no 



DD{^] -2Dr(^\+RR 



no 



(3) 



where DD, DR, and RR are pair counts of galaxies in 
the data-data, data-random, and random-random cata- 
logs with separation 9 + 69, and njj and Un are the mean 
number densities of galaxies in the data and random cat- 
alogs. This estimator has shown to be relatively insensi- 
tive to the size of the random catalog and handles edge 
corrections well (fKerscher. Szaoudi. fc Szal'avll2000() . We 
measure uj{9) in logarithmic bins of 0.2 in log(6') with 9 
measured in degrees. 

Measurements of uj{9) are known to be low by an ad- 
ditive factor known as the "integral constraint" (IC, see 
e.g. Peebles 1980) which results from using the data 
sample itself to estimate the mean galaxy density. This 
correction becomes important on scales comparable to 
the survey size, here equal to 0.5 degrees, the short axis 
of a single CFHT12k pointing. We therefore constrain 
our fits to a maximum 9 of 0.05 degrees, where the IC is 
negligible. 

5.2. Galaxy Clustering as a Function of Iab Magnitude 

Figures Eland [3 present uj{9) measured for independent 
magnitude bins in the range Iab — 18 — 24 on scales of 
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7" — 6'. The amplitude of the clustering is seen to 
smoothly decrease with magnitude, with fainter magni- 
tude bins showing less projected clustering. In Figure 
one can see that uj{6) approaches zero on large scales, 
while Figure shows clearly the effects of the IC; on 
scales larger than ~ 3' the amplitude of U!{9) begins to 
fall below the power-law form seen on smaller scales. In 
each of these plots we have measured the mean value of 
uj{0) in each bin. The error bars on uj{9) are calcu- 
lated from the standard error across our 15 independent 
pointings and should be a robust estimate of the error 
due to cosmic variance. We find that uj{0) measured by 
first summing the DD, DR, and RR counts over all the 
pointings, and then using Equation to estimate the 
clustering amplitude, is consistent with the mean values 
shown here. 

For each magnitude range we fit a power-law form, 
io{9)= Ai^9^, on scales ~ 7" — 3' and solve for the slope 
and amplitude using a least-squares fitting algorithm 
where we minimize x-squared, taking into account the 
covariance between bins as determined from bootstrap 
resampling. We first assume a slope of (5 = —0.8 to com- 
pare with other published results. TableEllists the ampli- 
tude of a;(6') measured at 1', A^{1'), in each magnitude 
bin for a slope of S — —0.8, along with an error from the 
least-squares fitting and the median number of galaxies 
within a single pointing in that magnitude range, equal 
to the median sample size. We find that the errors esti- 
mated from bootstrap resampling are nearly identical to 
those calculated above. In Figuredthese fixed-slope fits 
are plotted as dashed lines. As can be seen in this figure, 
a slope of i5 = —0.8 is generally consistent with the data 
for the magnitude ranges and scales we study here. 

In FigurelHlwe plot the amplitude of uj{9) at 1', A^{1'), 
as a function of median /ab magnitude for an assumed 
slope of (5 = -0.8. Our values agree well with the re- 
cent results of iMcCracken et all 1)20011) (who do not use 
independent magnitude bins, but always keep a bright 
limit of /ab = 18.5 and v ary the faint limit of their 
magnitude range ) and with 'WilsoJ ("2003*), though our 
data points are slightly lower (see Wilson (2003) for a 
figure comparing their data with other results in the lit- 
erature). However, our re s ults d o not agree well with 
IRraiTierd. Sma,i]. fc Monldl lllfltM. IWoods FahlmanI 
()1997|) . or iPostman et all l)1998f) . who all find signifi- 
cantly higher amplitu des at faint magnitudes, ^a b > 22 
The surveys o f both Brain grd. Small, fc Mo^ild l|1995D 
and IWoods fc Fahlman (199"7l) are aiiite small, covering- 
less than 150 arcmin^ each, and are therefo r e sub ject to 
significant cosmic variance. IPostman et al.l l|1998t) cover 
a very wide area (16 degrees^) separated over 256 con- 
tiguous individual pointings. We find that their values of 
10(9) are a factor of two higher than our measurements at 
/ab ^22. It is unlikely that this large of a discrepancy 
could be caused primarily by errors on their estimates of 
the IC, which would only affect larger scales than we are 
considering here, or by errors in stellar contamination 
estimates, which would have to be unreasonably large 
at these magnitudes. A much more likely explanation 
is that these differences may be the result of calibration 
errors in their photometry, as any zero-point differences 
across their pointings would lead to a systematic increase 
in the measured clustering on all scales, as pairs across 
pointings are used on even the smallest scales. This 



would also have a greater effect at fainter magnitudes, 
where the relative photometric error is larger. Here we 
measure ij;(0) only within individual pointings and do not 
attempt to compute lo{9) across pointings. The contin- 
ued decrease in clustering w hich we fi nd for fainter galax- 
ies is also seen by fMcCracken et al.l ((2001 ) and WilsonI 
(2003). 

As can be seen in Figure [7| the canonical slope of 
S ~ —0.8 does not the fit the data exactly at all mag- 
nitudes, and so we next relax the assumption of a fixed 
slope and measure both the slope and amplitude oiuj(9) 
for each magnitude bin on scales of ~ 7" — 3'. The result- 
ing slope as a function of median /ab magnitude is shown 
in Figure 1^1 and listed in Table |21 along with the ampli- 
tude, A^(l'). The slope varies between 6 ~ —0.7 1.0, 

with a maximum between /ab ~ 19 — 21 and shallower 
values at fainter magnitudes. As sh own in Figure El 
our values are consistent with those of IMcCracken et all 
(200l|, and the general trend we find in how the slope 
varies with magnitude is consistent with IPostman et alJ 
|l998), t hough th ey measure shallower slopes at all mag- 
nitudes. IWilsonI l)2003j) report that a slope of (5 = —0.8 
fits their data well for the magnitude range /ab = 20—24. 

A decrease in the slope of uj{9) at fainter magnitudes 
implies evolution in the slope of £,{r), such that the slope 
is shallower at higher redshift. Semi- analytic models pre- 
dict a decrease in the slope of f (r) from ~ 1.85 at z ~ 
to ~ 1.6 at 2; ~ 1 ijKauffmann et alJll999D . and measure- 
ments of ^(r) in early data from the DEEP2 G alaxy Red- 
shift Survey have a slope of 7 = —1.66 ±0.12 l|Coil et alJ 
^04). These results presented here are consistent with 
a similar decease in the slope at z ~ 1 (see Section 5.4 
for a discussion of the redshift distribution as a function 
of magnitude). 

5.3. Modelling uj{9) 

As the two-dimensional galaxy clustering seen in 
the plane of the sky is a projection of the three- 
dimensional clustering, ljj{9) is directly related to its 
three-dimensional analog ^ (r) . For a given ^ (r) , one can 
predict the amplitude and slope of uj{6) using Limber's 
equation, effectively integrating ^(r) along the redshift 
direction (e.g. Peebles 1980). If one assumes ^(r) (and 
thus Lo(9)) to be a power-law over the redshift range of 
interest, such that 



e(r, z) = 



then 



w{9) = 



r[(7-l)/2] A 



07- 



(4) 



(5) 



r(7/2) 

where F is the usual gamma function. Here we use the 
convention of r being a physical length, and x being a 
comoving length, such that r = ax, where a is the scale 
factor. We therefore quote the comoving scale-length of 
clustering as a function of redshift as xq {z) and the local 
value as a;o(0) (equal to ro(0)). The amplitude factor. A, 
is given by 



A 



(6) 



where dN/dz is the number of galaxies per unit redshift 
interval and g{z) depends on 7, and the cosmological 
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model: 

9iz) = (J) ^^'-'^Fix). (7) 

Here x is the comoving distance at a redshift z and F{x) 
is the curvature factor in the Robertson- Walker metric, 

ds^ ^ c^dt^ ^a^[dx^ lF{xf +x^de^ + x^sin^ed(j)% (8) 

If the redshift distribution of sources, dN/dz, is well- 
known, then the amplitude of can be predicted for 
a given power-law model of S,{x), such that measurements 
of to (9) can be used to place constraints on the evolution 
of^(x). 

5.4. Redshift Distribution as a Function of Magnitude 

In the literature there is a wide variety of adopted red- 
shift distributions for magnitude-limited samples, which 
is alarming as the predicted amplitude of uj{9) is quite 
sensitive to the assumed shape, and particularly the 
width, of the redshift distribution. Most authors use 
either galaxy evolution models or small observational 
samples (N^ 100) to construct dN/dz distributions. 
IMcCracken eTall l|200l use their own galaxy evolution 
models which are tested against an HDF sample of 120 
galaxies and which reproduce the trend of median mag- 
nitude with z and dispersion in redshift for a given mag- 
nitude range as seen by the Canada- France Redsh i ft Sur- 
vey (CFRS). iCabanac. de Laoparent. fc H icksorf (|200(]n 
use the CFRS luminosity function co mbine d with 
the P EGASE spectral atlas ijFioc fc Rocca-Volmerangd 
119971) to construct models of the dN/dz distribu- 
tions. As they point out, the choice of which at- 
las to use can affect the K-corrections by as much as 
50%, so that their d N/dz distribution is q uite sen- 
sitive to this choic e. Roche fc EalesI l|1999() use the 
iCharlot. Worthev. fc Bressan ( 1996) spectral evolution 
models combined with the luminosity function of the Las 
Campanas Redshift Survey data. 

Other authors have used observed redshift distri- 
butions from relatively small, deep spectroscopic sur- 
veys. Postman et al. ( 1998J use the empirical function 
dN/dz cx z^ex.p[—{z/zoy], fit to the CFRS spectroscopic 
data. These fits are presumably based on the sample of 
~ 600 galaxies for which CFRS obtained secure redshifts, 
with the highest z 1.3. They also use various evolving 
Schechter luminosity functions to model dN/ dz and find 
parame t ers wh ich provide good fits to their uj{9) results. 
iWilsonI (|200,1D use the form dN/dz oc z^exp(-z/zo), 
where the mean redshift is Szq and the median redshift 
is 2.67^0. They fit for zq using the SSA22 field sample 
of ICowieT Songaila. fc Bargcr (1999), which has ^ 200 
galaxies to / = 23.25, where corrections are made for 
faint galaxies without redshifts which are thought to lie 
between z = 1.5 — 2.0. This par ameterization of dN/d z 
is markedly different from that of'P ostman et alJ l)1998D , 
with more galaxies at higher redshift, and leads to very 
discrepant predictions of the amplitude of w(0) as a func- 
tion of magnitude (at /ab — 23 the CFRS dN/dz model 
predicts a value of A„ which is 6% lower for e = 0, 14% 
lower for e = 1, and 23% lower for e = 2; see Section 5.5 
for a definition and discussion of e) . 

Here we use a uniform sample of 2954 galaxies with 
secu re redshifts in early DEEP2 Galaxy Redshift Sur- 
vey ijDavis et al.ll2002l: iCoil et al.ll200l iMadewick et all 



12003ft data in the Extended Groth Strip (EGS), where 
redshifts are observed from z ~ — 1.5. We find that the 
form used bv IWilsonI l)2003l) . dN/dz cx z2exp(-z/zo), 
fits the DEEP2 data reasonably well. The overall nor- 
malization of dN/dz is irrelevant in estimating the am- 
plitude of oj{9)] it is the shape, and particularly the 
width, of dN/dz that is i mportant for ou r purp oses. The 
dN/ dz parameterization iPostman et alJ l)1998|) used for 
the CFRS data does not fit the DEEP2 sample; their 
narrower function contains too few objects at high red- 
shifts. 

To reconstruct the redshift distribution of sources in 
the DEEP2 data we have inversely weighted each galaxy 
by its p robability of selection , which is a function of mag- 
nitude f Newm an et al.ll2004(l . Figure 1101 shows the red- 
shift distribution of sources in the DEEP2 data along 
with our best-fit dN/dz distribut ion for diffe r ential mag- 
nitude bins. As discussed in iCoil et aTl l)2004f) . the 
DEEP 2 Galaxy Redshift Survey does not measure red- 
shifts for galaxies at z > 1.45 as the A 3727 A [Oil] 
doublet and Ca H-l-K absorption features move out of 
the observable spectral wavelength range. Therefore, 
any objects at z > 1.45 are not included in the sample 
used to derive dN/dz. This will affect the fainter magni- 
tude ranges in particular, which will be missing objects 
at z > 1.45. Chuck Steidcl has followed up our blue 
{{B — R) < 0.4) failures using LRIS-B (private commu- 
nication) and has found that of the 20 objects he has tar- 
geted, all were at z > 1.45, with a median z = 1.75 and 
a standard deviation of cr '--^ 0.30. This redshift distribu- 
tion is very similar to that shown injSteidel et al. (200j|) 
for the 'BM' galaxies, which lie at < z >= 1.70 ± 0.34, 
and it seems quite likely that our blue failures are in 
the 'BM' population. We correct our redshift distribu- 
tions by assigning redshifts to the blue galaxies without 
measured redshifts such that they have a Gaussian dis- 
tribution with a mean of z = 1.75 and a standard de- 
viation of a = 0.30. Our resulting redshift distributions 
are shown in Figure ITUl Fits for zq for each magnitude 
range are shown in Table 13 with and without the cor- 
rection for missed galaxies at z > 1.45, which affected 
only the fainter magnitude bins. Errors on zq were cal- 
culated using bootstrap resampling of the data and may 
underestimate errors due to cosmic variance, as the 2954 
redshifts which we use are from a single ~ 0.5 x 0.3 de- 
gree field. 

After fitting our dN/dz model to the DEEP2 spectro- 
scopic sample, we find that zq is well-fit as a linear func- 
tion of the median Jab magnitude: zq = —0.84 -I- 0.050 
median /ab? where the median redshift is 2.67zo and 
the mean redshift is 3zo. Figure 1111 plots our data 
points for the median redshift as a function of median 
/ab magnitude, with a solid line sh owing the linea r fit. 
We also include rece n t results from IWilsonI l)2003t) and 
iBrodwin et"aD l|2003D . IWilsonI (|2003D use the same func- 
tional form for dN/dz that we use here, though their fits 
of Zq vs magnitude do not match ours well at magnitudes 
fainter than /ab ~ 22, where they find a higher median 
redshift than we do. However, their results are based on 
a small sample of ~ 200 galaxies. Our results agree bet- 
ter with those of lBrodwin et al.l l)2003D . though they find 
systematically lower values for the median redshift as a 
function of magnitude and use a narrower dN/ dz fitting 
form which has fewer galaxies at higher redshifts. 
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We also include in Table 01 fits for dN/dz in integral 
magnitude ranges for both Rab and /ab bands, which 
are of interest for other analyses such as weak gravita- 
tional lensing studies. For the integral magnitude bins 
we fit for two distributions, dN/dz oc z^exp(— z/zq) and 
dN/dz cx 2:^exp(— z/zo)^'^, and find that they both fit the 
data well to within the noise. For the z^exp(— z/zq)^'^ 
distribution, the mean redshift is 2.09zo and the median 
1.91zo. Beyond our limiting magnitudes of -Rab = 24 
and /ab = 23.5, we extrapolate to fainter magnitudes 
using a simple linear fit to the zq values in the brighter 
bins. Note that all magnitudes are AB in the native 
CFHT 12k filter system. 

5.5. Constraining Evolution of Clustering 

Given the redshift distribution of sources, one can use 
measurements of u>{9) to either determine ^{x) at a given 
redshift, if the galaxies lie in a narrow redshift range, 
or parameterize the redshift-dependence of ^(a;, z) and 
solve for the local value of ^(x,0) and its evolution with 
r edshift. 

iGroth fc Peebles! ljl977|) suggested the 'e-model' for the 
evolution of the two-point correlation function: 

xo(z) = a;o(0)(l + z)-(3+^-^)/^ 



C(a;,z) 



(1 + ^) 



-(3+e-7) 



(9) 
(10) 



where 2; is a comoving distance, a:;o(0) is the local value 
of xo(z) and 7 is the slope of ^(x). In proper coordinates, 
the scale-length evolves as ro(z) = ro(0)(l + z)^''^^^'^^"' 
and ^(r) evolves as £,{r,z) = £.{r,Q){l + z)-(3+^). This 
model was motivated such that e = for clustering fixed 
in proper coordinates, where =constant and n is the 
mean galaxy density. Clustering which is fixed in comov- 
ing coordinates requires 6 = 7—8, which equals —1.2 for 
7 = 1.8 (equivalent io 5 — —0.8). Any value of e > 
leads to growth of clustering in proper coordinates. Lin- 
ear theory predicts e = 7 — 1, equal to 0.8 for 7 = 1.8, 
for growth of dark matter clustering in an = 1 cos- 
mology. Using this e parameterization for the growth of 
clustering, equation becomes 

/-.g(0)(l + z)-(3+^-^)g(z)(^)^dz 

Figure ^1 shows models of the predicted amplitude 
of Lo{0) as a function of magnitude for different values 
of a;o(0) and e, given our fits for dN/dz. An increase 
in xo(0) has the effect of shifting the predicted values 
higher, increasing the amplitude of clustering at all mag- 
nitudes, while increasing e creates a steeper slope in the 
evolution of clustering. We find that our data is not well- 
fit by any single value of e; the brightest galaxies are 
best-fit by e - and a;o(0) ~ 5 h'^ Mpc, to /ab ^ 21. 
The amplitude oiu){9) at 1' then begins to decrease more 
rapidly than an e = model predicts, such that the 
fainter galaxies imply a value of e ~ 3. This implies 
significant growth of clustering at intermediate redshifts 
(z ~ 0.5 — 1) and clustering fixed in physical coordinates 
in the more recent past (z ~ — 0.5). A model with 
e ~ 3 does not match our data within the errors across 
the full magnitude range; in particular, the model pre- 
dictions are low for /ab ~ 20 — 23. The e model is known 



not t o fit angular c lustering data well l|McCracken et all 
l2?)0TllWilsonl2nnl and is in fact not a physically realistic 
model given the presence of a galaxy bias which evolves 
with redshift. Therefore we attempt to find a new model 
which better fits the data. 

Given that uj{9) declines monotonically with the me- 
dian /ab magnitude of a sample, and that zq is well- 
fit as a linear function of the median /ab magnitude, 
we propose a linear model for the evolution of cco(z), 
2:0(2) = a^o(0)(l — Bz). We stress that this linear model 
for the evolution of a;o(z) is valid only in the redshift 
range this dataset is sensitive to, namely z ^ — 1.5, and 
a;o is expected to increase at hig her redshifts, as seen in 
both data and simulat ions (e.g.. lAdelberger et al.ll200'3 
iKauffmann et al1ll999() . For this model, ^(r) evolves as 

^(x,z)=[^1'^^N")(^-^-f =C(x,0)(l-i^2)- 
[ a; J [ X J 

(12) 

Replacing the xq{z) term in Equation with 2;o(0)(l — 
/3z), we can fit for B and xo(0). Fits for the linear model 
are shown in Figure El the best fit is xo(0) ^ 5 — 6 h^^ 
Mpc and B ~ 0.6 — 0.9. These values imply that a;o(z = 
0.8) - 1.4 - 3.1 h-'^ Mpc and a;o(z = 1.1) - 0.1 - 2.0 
Mpc. ^(r) results from early DEEP2 data measure 
xq{z = 0.8) = 3.5 ± 0.8 h'^ Mpc a nd xo(z = 1.1) = 
3.1 ± 0.7 /i-i Mpc ijCoil et al.ll2004|) . somewhat higher 
than what is implied here, though within 1-2— ct of the 
upper range of values. 

We do not necessarily expect the results to match ex- 
actly, as the ^(r) results are for a data sample with a lim- 
iting magnitude of Rab = 24.1, which is roughly half a 
magnitude brighter than our faintest uj{6) sample. Addi- 
tionally, a luminosity- dependent bias, which is expected 
to ex:ist at z ~ 1 fe.g.. lCoil et al.ll2'004HKauffmann et alJ 
ITSi), would cause the uj{9) results from our faintest 
magnitude bin to imply a somewhat lower value of a;o(z). 
There are other possible factors that may contribute 
to this discrepancy; a wider dN/dz distribution, espe- 
cially at the higher redshifts, would imply larger values 
of xo{z), though comparisons to other surveys do not 
indicate that our redshift distributions are too narrow. 
Limber's equation also has implicit assumptions which 
may break down at some level; for example. Limber's 
equation holds that the probability of finding a galaxy 
pair at a given separation depends on ^(r) multiplied 
by the selection function for each galaxy. This implies 
that galaxies have a selection function which is indepen- 
dent of their clustering. This assumption does not hold 
if clustering depends on environment or luminosity, for 
example, which is known to be true both locally and at 
intermediate redshifts. This may affect our ability to in- 
terpret Lu{9) results when using a sample that includes a 
mix of galaxy types, luminosities and redshifts. 

6. GALAXY CLUSTERING AS A FUNCTION OF {R - I) 
COLOR 

We now parallel the previous section, discussing the 
clustering as a function of color instead of magnitude. 
We plot in Figure E| the amplitude of uj{9) at 1' for 
nine independent {R — I) color bins, where we fit for 
the slope and amplitude simultaneously. Each bin has 
a width of (/? — /) = 0.2, and we use all galaxies in 
the magnitude range /ab = 18 — 24. Table [S] lists the 
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computed amplitude and errors of lu{9) measured at 1' 
for each color range. Similar to as a function of 

magnitude (Section 5.2), we find that the results are 
consistent if we first sum the DD, DR, and RR pairs 
across all the pointings before taking the ratio in Equa- 
tion (2) and errors are estimated using bootstrap resam- 
pling of the data. There is a strong, monotonic pro- 
gression between 0.2 < {R — I) < 1.5 of redder galax- 
ies having a higher clustering strength. This trend is 
not surprising, as red galaxies have long been known 
to be more strongly clustered than blue galaxies locally 
and this trend appears to continue to z ~ 1 (recent re- 
sults include Willmer et al. 1998, Zehavi et al. 2002, 
Coil et al. 2003). We also find that there is a sharp 
rise in the clustering amplitude for the bluest galaxies 
in our sample (—0.2 < {R — I) < 0.0), which have a 
clustering strength comp arable to the reddest g alaxies 
(1.2 < (R-I) < 1.6). ILandv. Szalav. fc Kod llT99l 
first reported such a marked correlation of clustering 
with color, in a sample of ~ 3000 galaxies selected by 
U — Rp color, where they found that the bluest and 
reddest galaxies showed much strong clustering than the 
sample as a whole. Mo re recently, with larger surveys, 
iMcCracken et al.l \2mi\\ and Wilson ( 200|) also find sim- 
ilar clustering trends with V — I color and report that 
the bluest and reddest samples are highly clustered. We 
note that the strong rise in clustering amplitude seen for 
the bluest galaxies is only apparent after rather stringent 
color cuts are made {R—I < 0.2), and that galaxies which 
are selected to be blue based on a _B magnitude alone do 
not show strong clustering (e.g., Efstathiou et al. 1991). 

Figure ^] shows the best-fit slope for each color bin, 
as listed in Table We find that galaxies with 0.2 < 
(R-I) < 1.0 have slopes of -0.6 < ^ < -0.8, while 
the bluest and reddest samples have steeper slopes, with 
-1-0 ^ <5 ^ -1.4. This implies that the slope of S,{r) for 
the bluer and redder samples is steeper than for galaxies 
with 0.2 < (i? - /) < 1.0. Locally, the slope of ^(r) in 
SDSS galaxies is found to be 7 = 1.86 ± 0.03 for redder 
galaxies with a restframe color of u* — r* > 1.8 and 

i ' = 1.41 zfc .04 for bluer galaxies with u* — r* < 1.8 
Zehavi et al. 2002 !). At interm ediate redshifts (0.1 < 
z < 0.5) [Shepherd et al.l 1)20011) find that galaxies with 
early-type SEDs have a clustering slope of 7 = 1.91±0.06 
while late-type galaxies show 7 = 1.59±0.08. Our results 
imply a slope of 7 > 2.0 for the reddest galaxies, and 
7 ~ 1.6 — 1.7 for typical blue galaxies. The very bluest 
galaxies, with {R — I) < 0, appear not to be typical 
intermediate-redshift blue galaxies, as described in the 
following sections. 

6.1. Redshift Distributions of Different Color Ranges 

Using the sample of 2954 early DEEP2 galaxies with 
secure redshifts between < 2: < 1.5 in the EGS we con- 
struct redshift distributions for sources as a function of 
{R — I) color as shown by the solid lines in Figure ^1 
Note that these histograms do not have the same scale 
on the y-axis, as the number of galaxies in each bin varies 
widely. The reddest galaxies, with 1.2 < (i? — /) < 1.6, 
lie in a relatively narrow redshift range centered near 
z ^ 0.85. Successively bluer bins have both a lower me- 
dian redshift and a wider redshift distribution. Galaxies 
with 0.0 < (R — I) < 0.4 appear to be predominately 
centered at z ~ 0.35, though there is a tail to higher 



redshifts. The bluest galaxies (-0.2 < (i? - /) < 0.0) 
have a narrow redshift distribution between z ~ — 0.5. 
As discussed above, the DEEP2 Redshift Survey can not 
measure redshifts for galaxies at z > 1.45 as observable 
emission and absorption line features move out of the 
spectral wavelength range of our data. These redshift his- 
tograms are therefore missing galaxies at z > 1.45. The 
redshift success rate of the DEEP2 survey in the EGS is 
fairly insensitive to color in the range 0.2 < (R—I) < 1.6, 
where the success rate is > 80% and where 90% of the 
targeted objects lie, but it begins to drop blueward of 
{R-I) 0.2. For galaxies with 0.0 <{R-I) < 0.2 the 
redshift success rate is ~ 55%, while the bluest galax- 
ies, with —0.2 < {R—I) < 0.0, have a success rate of 
^ 30%. We believe that most of the observed galaxies 
without redshifts lie at z > 1.45, as these blue galax- 
ies should have emission lines which we would detect at 
lower redshifts. The sample of 20 DEEP2 blue galaxies 
without measured redshifts which Chuck Steidel has fol- 
lowed were all found to lie at z > 1.45 (see Section 5.4 
for details). We therefore correct the bluest two redshift 
distributions by assigning redshifts to the galaxies with- 
out measured redshifts such that they have a Gaussian 
distribution with a mean of z = 1.75 and a standard de- 
viation of cr = 0.30, matching the distribution of sources 
followed by Steidel. These corrected redshift distribu- 
tions can been seen in Figure IT?^ as the dashed lines for 
the two bluest distributions. 

6.2. Clustering Scale-length as a Function of {R — I) 

Color 

Figure lTTI shows the comoving scale-length xq as a func- 
tion of {R—I) color calculated using the observed redshift 
distributions shown in Figure ITHI in solid lines, assuming 
that Xq is constant over the redshift range of galaxies in 
each color bin. The errors on xq are derived from the 
errors on lu{9) alone and do not include uncertainties in 
the redshift distribution; therefore they may be under- 
estimates of the true error. Values of xq and associated 
errors are given in Table |5| As dividing the galaxies into 
color bins results in more uniform populations within a 
single bin, we believe that Limber's equation should be 
more robust in interpreting w(0) as a function of color. 

For the reddest galaxies, with 1.4 < {R — I) < 1.6, 
wc find xo — 5.02 ± 0.26 Mpc (note that this is 
the comoving scale- length at the redshift of the sample) . 
These galaxies, centered at z ~ 0.85, have absorption- 
dominated spectra and are likely to be massive pro- 
genitors of local elliptic als. Our value is somewhat 
lower than that found bv lBrown et alJ l)2003j) . who mea- 
sure early-type galaxies with Bw — R > 1-28 between 
0.7 < z < 0.9 t o have xo = 6.6 ± 0.8 h'^ Mpc. Locally, 
iBudavarT eTal] l|2003|) finds from an angular clustering 
study of SDSS galaxies that a;o(0) = 6.59 ± 0.17 /i"! 
Mpc and 5 = —0.96 ± 0.05 for early-type galaxies with 
—21 > Mr > —23. Our results at z '--^ 0.85 are con- 
sistent with predictions of little evol ution of xn(z) for 
early-type galaxies from z ~ 1 to ifKauffmann et afl 
IT999tlBenson et al.ll200l . Our reddest sample may be a 
less-extreme version of the extremely red object (ERO) 
population which are seen to have large clustering at 
z ^ 1 with xn 7.5 - 12 /i-i Mpc (e.g., IDaddi et alJ 
l200ltlFirth etal]l200l . 

The comoving scale-length decreases as galaxies be- 
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come bluer to a minimum of xq — 1.60 ± 0.23 h^^ Mpc 
at 0.2 < (i? — /) < 0.4. These galaxies lie predomi- 
nately between z ~ 0.3 — 0.6, though there is a small 
tail at z > 1.0. The redshift distribution of galaxies with 
0.0 < {R — I) < 0.2 is strongly bimodal; this sample 
appears to be a mix of two populations: relatively local 
galaxies similar to the sample with 0.2 < (i? — /) < 0.4 
and a z > 1.4 population of bright, star-forming galax- 
ies. If the angular clustering is due entirely to the lower- 
redshift population, then xq = 2.47 ± 0.32 h^^ Mpc, 
where we have used the uncorrected redshift distribution 
shown as a solid line in Figure^] If the clustering is due 
entirely to the higher-redshift population, then xq = 3.8 
Mpc. More likely, however, is that the observed uj{0) 
value is a mixture of these two populations, and we would 
like to separate their contributions to ll!{9) and measure 
values of xq for each population, as they may have differ- 
ent intrinsic clustering properties. It can be shown that 
to (9) from a combined sample of disparate populations 
with different clustering properties and non-overlapping 
dN/ dz distributions adds linearly such that 

C0{0)total = fluJi{B) + (1 ~ hfu2{e) (13) 

where /i is the fraction of the total number of galax- 
ies in one population, (1 — /i) is the fraction in the 
other population, and uii(9) and uj2{9) are the an- 
gular clustering amplitudes of each population (e.g., 
iMeiksin. White, fc Peacockll999() . The value of xq = 3.8 
Mpc which we find for the z > 1.4 population is 
therefore a lower limit on the clustering scale-length of 
that population, as some of the observed projected clus- 
tering is due to the lower-redshift sample. If we use the 
value of uj{9) observed for the 0.2 < {R — I) < 0.4 sam- 
ple as the clustering contribution from the lower-redshift 
population, we can then subtract this from the observed 
u){9) amplitude for the 0.0 < [R — I) < 0.2 sample, as- 
suming that ~50%-75% of the galaxies in this color bin 
are in the lower-redshift population. Using the redshift 
distribution of the z > 1.4 galaxies shown in Figure [TBI 
we find xq = 7.8 ± 0.3 Mpc for the z > 1.4 galaxies 
if 50% of the galaxies in this color bin are in this pop- 
ulation (corresponding to all the galaxies without mea- 
sured redshifts), and xa = 5.2 ± 0.4 h'^ Mpc if 25% 
of the galaxies are at z > 1.4 (corresponding to half of 
the galaxies without measured redshifts). These values 
most likely bracket the true value, as the percentage of 
galaxies without measured redshifts at z > 1.4 is likely 
to be between 50%-100%. These results are shown in 
Figure El It is worth noting that this error does not in- 
clude uncertainties in the redshift distribution of z > 1.4 
galaxies. 

In the bluest color bin, -0.2 < {R - I) < 0.0, the 
number of galaxies drops dramatically and the redshift 
distribution for objects at z < 1.4 appears to be relatively 
narrow, though the robustness of this width is uncertain 
given the small sample size. If the clustering of the bluest 
objects is due to galaxies with the uncorrected redshift 
distribution shown as a solid line in the upper left panel 
of Figure [TEl then xq = 3.69 ± 0.87 h-^ Mpc. However, 
as the redshift success rate is only ~ 30% for galaxies in 
this bin in the DEEP2 EGS region, it is likely that most 
of the objects in this color range are at z > 1.4, since they 
must be actively forming stars to have such blue colors 
and would therefore have strong nebular emission lines 



which would allow DEEP2 to determine their redshifts 
easily if at z < 1.4. A third population also contributes 
to objects found in this color bin: broad-line AGN found 
in the DEEP2 sample disproportionally have (R—I) < 0. 
In the EGS, where the photometric-redshift pre-selection 
of targets is not apphed, 11% of the galaxies in the bluest 
bin are identified as broad-line AGN, whereas < 1% are 
objects with {R — I) > are so classified. These per- 
centages are lower limits, since the spectroscopic signa- 
tures of weak AGN could be dominated by galaxy light 
and therefore remain undetected. The AGN identified in 
DEEP 2 spectra span the full redshift range < z < 1.5. 
The bluest color bin, -0.2 < {R- 1) < 0.0, thus includes 
at least three distinct populations: faint blue galaxies at 
low redshift, bright star-forming galaxies at z > 1.4, and 
broad-line AGN with strong, blue continua. 

Separating the contributions to the observed total ol!{9) 
from each population is not entirely straightforward, as 
there are uncertainties in both the estimates of the frac- 
tion of objects in each population and their redshift dis- 
tributions. However, given the small percentage of local 
galaxies seen in the observed redshifts for this bin and 
the low clustering amplitude of the local star-forming 
galaxies in the 0.2 < {R — I) < 0.4 bin, which have 
Xq = 1.60 ± 0.23 h'^ Mpc, it seems unlikely that local 
galaxies alone account for the large clustering signal seen 
here. A sample of local galaxies with the same relatively 
redshift distribution and clustering scale-length as those 
in the 0.2 < {R — T) < 0.4 color range, would have a 
uj{9) value 13% as large as what is measured in th is bin, 
in opp osition to the concl usions of^Land_y. Szalav. fc Kool 
l)1996l) . lMcCracken et all pOOl|) amf Mlso n (2003), who 
suggest that the strong clustering of blue galaxies is dom- 
inated by a local galaxy population. 

If the observed uj{9) in this color bin was due en- 
tirely to z > 1.4 galaxies, the clustering scale- length 
for that population would be xo — 7.5 ±1.5 h^^ Mpc. 
If AGN were the sole contributing population and had 
the same redshift distribution as the local galaxies in 
this color bin, they would have a clustering scale-length 
of Xo = 3.69 ± 0.87 Mpc, the same as for local 
galaxies. Observations of AGN and QSO clustering have 
found scale-lengths comparable to local early-type galax- 
ies, with Xq 6 — 8 Mpc, an d very little evolu- 
tion i n the scale-leng th with redshift ijCroom et alJl2003 
Grazi aiTet al.ll20M) . It therefore seems quite feasible 
that AGN contribute significantly to the large cluster- 
ing seen in this color range. Assuming that the observed 
Lo{9) is due to a combination of these three populations, 
and that each of them is totally independent, would lead 
to much larger values for xq of each subsample, possibly 
in excess of other estimates of the clustering of z > 1.4 
galaxies and AGN. It is likely that not only the autocor- 
relation of each population, but also its cross-correlation 
with the other sets of extremely blue objects, leads to 
the strong observed clustering in this color bin. 

In order to check the degree of independence between 
the redshift distributions in each sample, we measure the 
cross-correlation between galaxies in color bins i and j, 
Wij{9), and normalize each correlation amplitude, A^.., 
by the square root of the auto-correlation function am- 



12 



plitude, At^..: 

VA^„(1')A^„(1') 

where we have measured w{l') assuming a constant slope 
of 5 = —0.8 and Wij is defined by equation O where DD 
is replaced with D^EXi. This results in a covariance ma- 
trix, shown in Table |d1 that provides information on the 
physical covariance between galaxies in each color bin. 
Our results are generally consistent with the covariance 
of the redshift distributions seen in Figure E| Each 
color sample shows significant covariance with the im- 
mediately adjacent samples, and galaxies in the reddest 
four color bins, with 0.8 < {R — I) < 1.6, all have sig- 
nificant covariance with each other. The bluest color bin 
does not have significant covariance with galaxies with 
0.2 < {R — I) < 0.4, implying that the relatively local, 
blue, star-forming galaxies with low correlation seen in 
the 0.2 < {R — I) < 0.4 range contribute very little to 
the clustering seen in the bluest population, as we as- 
sumed above. However, the z > 1.4 galaxy population 
must contribute to both of the bluest two color bins. This 
covariance matrix provides a reassuring check on our red- 
shift distributions used to interpret uj{9) as a function of 
color. 

7. SUMMARY AND CONCLUSIONS 

We measure the projected angular correlation func- 
tion, uj{9), using deep photometric imaging with a sam- 
ple of 350,000 galaxies to Iab = 24, a significant sub- 
set of which is used as photometric data for the DEEP2 
Galaxy Redshift Survey. Our data covers a total of 5 
degrees^ over 5 separate fields, a larger area than has 
been studied by previous deep angular clustering sur- 
veys. This increase in area, combined with galaxy red- 
shift distributions derived from a relatively large, uni- 
form spectroscopic sample of 2954 DEEP2 galaxies be- 
tween < z < 2, leads to tighter constraints on models 
than has been previously available. Using robust redshift 
distributions as a function of color also allows for quan- 
titative interpretations of the strong trends seen in uj{0) 
vs. observed (i? — /). 

We find that there has been significant growth of clus- 
tering from z > 1 to z ~ 0, and that the 'e-model' pro- 



posed bv lGroth fc Feeble^ l)1977f) does not describe the 
evolution of clustering well. We propose an alternative 
model in which the comoving scale-length, xq{z), evolves 
linearly over the redshift range we are sensitive to here 
and find that this model better describes the data. We 
find that xq 1 — 3 h^^ Mpc at z ~ 1, though these re- 
sults are not straightforward to interpret, as they result 
not only from intrinsic evolution, but also any changes 
in the mean luminosity, color, and type of galaxies with 
redshift. 

We also find that uj{d) is a very sensitive function 
of observed {R — I) color, with both the reddest and 
bluest galaxies showing strong clustering amplitudes and 
steeper clustering slopes. Red galaxies with 1.4 < {R — 
I) < 1.6 have a comoving scale-length of xq = 5.02±0.26 
h^^ Mpc and lie in a narrow redshift range centered at 
z ~ 0.85. Their spectra indicate that these galaxies have 
older stellar populations and are likely to be progenitors 
of local ellipticals. The clustering amplitude decreases 
with color to a minimum at 0.2 < {R — I) < 0.4, where 
Xq = 1.60 ± 0.23 h^^ Mpc for blue star-forming galax- 
ies between z ~ 0.3 — 0.6. The bluest galaxies, with 
{R — I) < 0, are very highly clustered, and include addi- 
tional populations of z > 1.4 bright, star- forming galax- 
ies with a scale-length of xq > 5 Mpc and broad-line 
AGN. We find that the clustering of the bluest galaxies 
is not likely to be dominated by a local population or 
due entirely to high-redshift galaxies. 
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Table 1. Details of our CFHT12k observations. 



Field pointing RA Dec. exp. time (hrs) seeing date of obs. 

B R I (FWHM) 



0230 1 02 27 30 

2 02 30 00 

3 02 32 30 
1052 1 10 52 43 
1416 1 14 15 29 

2 14 17 55 

3 14 20 21 
1651 1 16 48 00 

2 16 51 30 

4 16 55 28 
2330 1 23 26 56 

2 23 30 00 

3 23 33 03 

4 23 25 40 

5 23 34 00 
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Table 2. Power law fits to uj{9)= A^9^ as a function of 7ab magnitude. 



^ab 


range 


median 7ab 


# gals/pt. 


{S = -0.8) 


Aw(l') 


s 


18 


- 19 


18.6 


152 


0.187 ± 0.023 


0.184 ±0.027 


-0.82 ±0.10 


19 


- 20 


19.6 


453 


0.070 ± 0.010 


0.067 ±0.010 


-1.01 ±0.09 


20 


- 21 


20.6 


1224 


0.060 ± 0.004 


0.045 ± 0.005 


-0.99 ± 0.05 


21 


- 22 


21.6 


2812 


0.031 ± 0.001 


0.027 ± 0.002 


-0.90 ± 0.04 


22 


- 23 


22.6 


5921 


0.015 ± 0.001 


0.017 ±0.001 


-0.74 ± 0.05 


23 


- 24 


23.6 


11879 


0.007 ± 0.001 


0.008 ± 0.001 


-0.70 ± 0.07 



Table 3. Fits for the redsiiift distribution of DEEP2 spectroscopic 
sources as a function of magnitude, for differential magnitude bins, 
assuming dN/dz oc z'^ exp{— z / zq) . The corrected values of zq include 

GALAXIES AT Z > 1.4 (SEE TEXT FOR DETAILS). 



/ab range uncorrected 20 corrected 20 zo error 



18 - 19 

19 - 20 

20 - 21 

21 - 22 

22 - 23 

23 - 24 



0.091 
0.136 
0.197 
0.237 
0.278 
0.285 



0.091 
0.136 
0.197 
0.239 
0.288 
0.324 



0.010 
0.007 
0.005 
0.005 
0.006 
0.010 



Neuschaefer, L. W., & Windhorst, R. A. 1995, ApJ, 439, 14 
Newman, J., et al. 2004, in preparation 
Norberg, R, et al. 2001, MNRAS, 328, 64 
Peebles, P. J. E. 1975, ApJ, 196, 647 

Peebles, P. J. E. 1980. The Large-Scale Structure of the Universe, 

Princeton, N.J., Princeton Univ. Press 
Postman, M.. Lauer, T. R., Szapudi, I., & Oegerle, W. 1998, ApJ, 

506, 33 

Roche, N., & Eales, S. A. 1999, MNRAS, 307, 703 
Roche, N., Shanks, T., Metcalfe, N., & Fong, R. 1993, MNRAS, 
263, 360 



Shepherd, C. W., et al. 2001, ApJ, 560, 72 

Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 
D. K., Reddy, N. A., & Hunt, M. P. 2004, ApJ, 604, 534 

Willmer, C. N. A., da Costa, L. N., & Pellegrini, P. S. 1998, AJ, 
115, 869 

Wilson, G. 2003, ApJ, 585, 191-204 
Woods, D., & Fahlman, G. G. 1997, ApJ, 490, 11 
York, D. G., et al. 2000, AJ, 120, 1579 
Zehavi, I., et al. 2002, ApJ, 571, 172 



14 



Table 4. Fits for the redshift distribution of DEEP2 spectroscopic 

SOURCES as a function OF MAGNITUDE, FOR INTEGRAL MAGNITUDE BINS. 



Rab range 






7ab range 




zo"" 


18 - 20 


0.104 


0.145 


18 - 20 


0.123 


0.171 


18 - 21 


0.136 


0.188 


18 - 21 


0.168 


0.232 


18 - 22 


0.168 


0.232 


18 - 22 


0.206 


0.284 


18 - 23 


0.207 


0.285 


18 - 23 


0.242 


0.332 


18 - 24 


0.255 


0.350 


18 - 24"= 


0.283 


0.389 


18 - 25^^ 


0.286 


0.392 


18 - 25^^ 


0.322 


0.443 


18 - 26"= 


0.324 


0.443 


18 - 26"= 


0.362 


0.496 


18 - 27" 


0.361 


0.494 


18 - 27^^ 


0.401 


0.550 



''for dN/dz oc z^exp(— z/zq) 
''for dN/dz cx z^exp(-2:/2o)^-^ 

^extrapolated using a linear fit (see text for details) 



Table 5. Power law fits to uj{6)= Ai^d^ as a function of (R — I) 
color, using the observed redshift histograms shown in figure 16 
(solid lines). 



(R- 


I) 


median 


# gals/pt 


Au,{l') 


Au,{l') S 


2:0(2) 


range 


(R-I) 




{S = -0.8) 




Mpc 


-0.2 - 


- 0.0 


-0.05 


274 


0.118 ± 0.021 


0.113 ±0.023 -1.36 ±0.21 


3.69 ±0.87 


0.0 - 


0.2 


0.14 


1869 


0.037 ± 0.003 


0.041 ±0.005 -0.69 ±0.11 


2.47 ± 0.32 


0.2 - 


0.4 


0.31 


5796 


0.012 ± 0.001 


0.014 ±0.001 -0.59 ±0.08 


1.60 ±0.23 


0.4 - 


0.6 


0.50 


5745 


0.019 ± 0.001 


0.021 ±0.001 -0.69 ±0.06 


2.78 ±0.36 


0.6 - 


0.8 


0.69 


4415 


0.031 ± 0.001 


0.032 ± 0.002 -0.78 ± 0.05 


3.34 ±0.30 


0.8 - 


1.0 


0.88 


2364 


0.048 ± 0.003 


0.049 ±0.003 -0.76 ±0.07 


3.38 ±0.22 


1.0 - 


1.2 


1.09 


1221 


0.080 ± 0.005 


0.072 ±0.006 -1.02 ±0.07 


4.18 ±0.22 


1.2 - 


1.4 


1.29 


861 


0.098 ± 0.006 


0.093 ±0.006 -1.25 ±0.06 


4.25 ± 0.33 


1.4- 


1.6 


1.46 


308 


0.171 ±0.014 


0.144 ±0.016 -1.14 ±0.10 


5.02 ± 0.26 



Table 6. Covariance of u!(9) as a fuxction of (R — I) color 



(R-I) -0.2-0.0 0.0-0.2 0.2-0.4 0.4^0.6 0.6-0.8 0.8-1.0 1.0-1.2 1.2-1.4 1.4^1.6 

range 



-0.2-0.0 


1.00 


0.66 


0.14 


-0.06 


-0.22 


-0.19 


-0.06 


-0.14 


-0.36 


0.0-0.2 


0.66 


1.00 


0.62 


0.12 


-0.13 


-0.16 


-0.19 


-0.20 


-0.26 


0.2-0.4 


0.14 


0.62 


1.00 


0.64 


0.26 


0.09 


-0.06 


-0.01 


-0.02 


0.4-0.6 


-0.06 


0.12 


0.64 


1.00 


0.70 


0.38 


0.22 


0.25 


0.10 


0.6-0.8 


-0.22 


-0.13 


0.26 


0.70 


1.00 


0.72 


0.57 


0.63 


0.43 


0.8-1.0 


-0.19 


-0.16 


0.09 


0.38 


0.72 


1.00 


0.87 


0.71 


0.49 


1.0-1.2 


-0.06 


-0.19 


-0.06 


0.22 


0.57 


0.87 


1.00 


0.92 


0.54 


1.2-1.4 


-0.14 


-0.20 


-0.01 


0.25 


0.63 


0.71 


0.92 


1.00 


0.81 


1.4-1.6 


-0.36 


-0.26 


-0.02 


0.10 


0.43 


0.49 


0.54 


0.81 


1.00 
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dx 

Fig. 1. — Residuals in the astrometric solution, with units in pixels. The rms is ~0.025 pixels, corresponding to ^ 0.005". 
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Fig. 2. — Differences in magnitudes for stars detected on multiple exposures. The rms at bright magnitudes is ~0.02. 
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Fig. 3. — Color-color distribution for compact galaxies and stars in one pointing; these distributions are used to define the probability 
that a given object is a galaxy. 




Fig. 4. — Probability distribution that a compact object is a galaxy for data in all pointings. The distribution is quite bi-modal; most 
compact objects are clearly classified as stars or galaxies. 
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Fig. 5. — Galaxy counts for -Bab, Rab, and Iab bands. The dotted lines are linear fits discussed in the text (see Section 4). 
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Fig. 6. — The angular correlation function, U!{0), measured for independent /ab magnitude bins, shown on scales ~ 7" — 6'. 
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Fig. 7. — The angular correlation function, i^^iff), plotted on a log scale for different /ab magnitude ranges. Power law fits for ui{8)- 
on scales ~ 7" — 2' are shown as dashed lines, assuming a constant slope of 5 = —0.8. 
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Fig . 8. — The amplitude of uj(d) measured at 1 ' shown as a function of the median /ab magnitude, assuming a slope of 5 = —0.8. Results 
from lMcCracken et al.l J2002) and lWilsonj Oooll are plotted as well, with dotted error bars for clarity. 
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Fig. 9. — The best-fit slop e, 5, of u!(d) as a func tion of median /ab magnitude. The slopes are fit on scales of ~ 7" — 3'. Results from 
IMcCracken et al.l <200ll) and lPostman et al.l 11993) are shown for comparison. 
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Fig. 10. — Redshift histograms of galaxies in the early DEEP2 Redshift Survey data for independent Iab magnitude ranges. Fits for 
dN/dz oc exp{— z / zq) are plotted as solid curves, and values for zq are listed in Table The mean redshift in each magnitude range is 
3^0 a-iid the median redshift is 2.67zo- 
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Fig. 11. — Median redshift as a function of /ab magnitude for data from the DEEP2 Redshift Survey, using the fits shown in Figure IT?)1 
where the median redshift is 2.67^o- The solid line is a linear fit to the data; zo = —0.84 + 0.050 median /ab- Results from two recent 
surveys are also shown. 
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Fig. 12. — The amplitude of uj{8) at 1' as a function of the median /ab magnitude compared with predictions for different values of the 
local clustering scale-length, xo{0), and an evolutionary parameter, e, such that xo(z) = xo{0){l + 2)"'^+'^ , 
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Fig. 14. — T he amplitude of lij{0) measured at 1' shown as a function of (R — I) color, where the slope has been measured simultaneously 
(see Figure [T^ . 
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Fig. 15. — The best-fit slope, <5, of i^{d) as a function of {R — I) color. 
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Fig. 16. — Redshift distribution of DEEP2 EGS sources as a function of (R — I) color. The daslied lines for the two upper left plots show 
the assumed redshift distribution for galaxies at z > 1.45 (see text for details). 
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Fig. 17. — The comoving scale-length 3:0(2:) as a function of observed (R — I) color. In addition to the scale-length derived using the 
observed redshift distributions shown in Figurc lTCl wc also plot values of the scale-length derived for 2: ~ 1.7 galaxies (see text for details). 



